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Abstract. We consider a class of supersymmetric inflation models, in which min- 
imal gauged F-term hybrid inflation is coupled renormalisably to the minimal super- 
symmetric standard model (MSSM), with no extra ingredients; we call this class the 
"minimal hybrid inflationary supersymmetric standard model" (MHISSM). The sin- 
glet inflaton couples to the Higgs as well as the waterfall fields, supplying the Higgs 
/i-term. We show how such models can exit inflation to a vacuum characterised by 
large Higgs vevs, whose vacuum energy is controlled by supersymmetry-breaking. The 
true ground state is reached after an intervening period of thermal inflation along the 
Higgs flat direction, which has important consequences for the cosmology of the F-term 
inflation scenario. The scalar spectral index is reduced, with a value of approximately 
0.976 in the case where the inflaton potential is dominated by the 1-loop radiative 
corrections. The reheat temperature following thermal inflation is about 10^ GeV, 
which solves the gravitino overclosure problem. A Higgs condensate reduces the cos- 
mic string mass per unit length, rendering it compatible with the Cosmic Microwave 
Background constraints without tuning the inflaton coupling. With the minimal U(l)' 
gauge symmetry in the inflation sector, where one of the waterfall fields generates 
a right-handed neutrino mass, we investigate the Higgs thermal inflation scenario in 
three popular supersymmetry-breaking schemes: AMSB, GMSB and the CMSSM, fo- 
cusing on the implications for the gravitino bound. In AMSB enough gravitinos can 
be produced to account for the observed dark matter abundance through decays into 
neutralinos. In GMSB we find an upper bound on the gravitino mass of about a TeV, 
while in the CMSSM the thermally generated gravitinos are sub-dominant. When Big 
Bang Nucleosynthesis constraints are taken into account, the unstable gravitinos of 
AMSB and the CMSSM must have a mass 0(10) TeV or greater, while in GMSB we 
find an upper bound on the gravitino mass of 0(1) TeV. 
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1 Introduction 

Inflation is the accepted paradigm for the very early universe, thanks to its power to 
account accurately for cosmological data in one simple framework. However, it raises 
a number of theoretical problems, principally the identity of the inflaton, the flatness 
of its potential, and how it is coupled to the Standard Model. 

A technically natural way of achieving a flat potential is through supersymmetry 
(SUSY). However, the flatness is generically spoiled in supergravity [1], which must 
be taken into account if the inflaton changes by an amount of order the Planck scale 
or more ( "large- fleld" inflation). Given the large parameter space of supergravity 
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theories, this motivates starting the search for a supersymmetric theory of inflation 
with small-field inflation, in the context of a renormalisable theory. 

At the same time, low energy supersymmetry remains an attractive theoretical 
framework in which to understand the smallness of the electroweak scale relative to 
the Planck scale. The Minimal Supersymmetric Standard Model (MSSM) is the most 
economical possibility to combine low energy SUSY with the phenomenological triumph 
of the Standard Model (although the high Higgs mass and the absence of positive results 
from the Tevatron and LHC increases the amount of parameter tuning required). 

Indeed, the MSSM itself can realise inflation along one of the many flat directions 
[2] with the addition of non-renormalisable couplings. Inflation takes place near an in- 
flection point in the potential, where trilinear and soft mass terms are balanced against 
each other, although the amount of tuning required [3] reduces the attractiveness of 
the scenario. The tuning can be reduced by extending the MSSM [4, 5]. 

The simplest class of renormalisable supersymmetric inflation models is minimal 
F-term hybrid inflation, by which we mean the first supersymmetric model of Ref. [1], 
characterised by the superpotential 

Wi = Xi(!?'^S - A'PS. (1.1) 

General theoretical considerations of small-field infiation drive one towards this model 
[6], which works without a Planck-scale infiaton field, non-renormalisable operators, 
or supersymmetry-breaking terms. It invokes an infiaton sector of (at least) 3 chiral 
superfields, consisting of the infiaton itself, S, and two waterfall [7] fields, $, $, with 
an optional gauge superfield.^ Because the the infiaton field appears linearly in the 
superpotential, it does not suffer from the generic supergravity problem of Hubble-scale 
mass terms during infiation [1]. 

In its standard form, however, F-term hybrid infiation suffers from a number 
of problems which reduce its power to fit cosmological data. First and foremost is 
the gravitino problem, which limits the reheat temperature to be unnaturally small 
compared with the infiation scale. Of less severity is the spectral index problem. If 
the infiaton potential is dominated by the 1-loop radiative corrections, F-term hybrid 
infiation predicts that the spectral index of cosmological perturbations e-foldings 
before the end of infiation is = 1 — 1/A^. For the canonical 60 e-foldings, this is more 
than Icr above the WMAP7 value Ug = 0.963 ± 0.012. Finally, many models generate 
cosmic strings, and the CMB constraints on their mass per unit length forces one to 
very weak infiaton couplings, where — 1 [8]. 

There also remains the question of how the infiaton sector is coupled to the MSSM. 
If we restrict ourselves to renormalisable theories combining minimal U(l)'-gauged F- 
term hybrid infiation with the MSSM, with no other fields, and preserving all the 
symmetries, the choices are limited. The singlet infiaton 5* can couple in the super- 
potential only to the product of the Higgs fields or the square of the right-handed 
neutrino fields (which we take to be included the MSSM). If the MSSM fields have 

"'^The number of chiral superfields can be reduced to 2 without a gauge field, or if they are in a 
real representation. 
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non-trivial charge assignments under the U(l)' of F-term inflation, the coupling of S 
to the neutrinos is forbidden, and its place taken by one of the waterfall fields. This 
has the nice feature of generating a see-saw mechanism, with the neutrino masses also 
controlled by the vev of the waterfall fields. Neutrino masses are also allowed if the 
waterfall fields are U(2) triplets, with SU(2)/j as a subgroup. 

We will refer to the minimal case where the symmetry of the waterfall fields is 
U(l)' as the Minimal Hybrid Inflationary Supersymmetric Standard Model (MHISSM). 
In the model, it is very natural that the gauge singlet infiaton S should be coupled 
both to the waterfall fields and to the Higgs fields, which mixes the standard MSSM 
Higgs flat direction with the hybrid inflation waterfall direction. If the coupling of 
the infiaton to the Higgs is smaller than to the waterfall fields, inflation ends with the 
development of vevs for the Higgs multiplets, /11.2, breaking the electroweak symmetry. 
Soft terms lift the flat direction, and if certain constraints are satisfied, the Higgs fields 
will finally reach the standard vacuum after a period of thermal inflation, with a reheat 
temperature of about 10^ GeV. This solves the gravitino overclosure problem, and Big 
Bang Nucleosynthesis constraints can be satisfied with massive (0(10) TeV or more) 
or stable gravitinos [9-12]. 

We call this second period of accelerated expansion Higgs thermal inflation. It 
is a natural consequence of the coupling of the F-term hybrid infiaton to the Higgs 
fields, and offers a generic solution to the gravitino problem. At the same time, a TeV- 
scale vacuum expectation value for the infiaton generates an effective /^-term. The 
model was first introduced in Ref. [13] in the context of Anomaly-Mediated Super- 
symmetry Breaking (AMSB). We termed the version of AMSB there deployed strictly 
anomaly mediated supersymmetry breaking (sAMSB), because D-terms associated with 
the U(l)' symmetry resolve the AMSB tachyonic slepton problem, without requiring 
an additional explicit source of supersymmetry breaking. 

In this paper we demonstrate that the interesting cosmological consequences, in 
particular Higgs thermal inflation, are a result of the structure of the model at the 
inflation scale, and not of the particular supersymmetry-breaking scenario. We derive 
the effective potential for the combination of fields driving thermal inflation, and the 
constraints on the soft breaking parameters for a phenomenologically acceptable ground 
state, in three popular supersymmetry-breaking scenarios: anomaly-mediated (AMSB), 
gauge-mediated (GMSB) and the constrained minimal supersymmetric standard model 
(CMSSM). We find that the lower reheat temperature following thermal inflation solves 
the gravitino problem in the CMSSM, while in AMSB enough gravitinos can be produced 
to account for the observed dark matter abundance through decays into neutralinos. 
In GMSB we find an upper bound on the gravitino mass of about a TeV, derived from 
constraints on NLSP decays during and after Big Bang Nucleosynthesis (BBN). 

F-term models with Higgs thermal inflation have other important features. The 
spectral index of scalar Cosmic Microwave Background fluctuations Ug is reduced, as 
fewer e-foldings of F-term inflation are required. In the range of couplings for which 
the 1-loop radiative corrections dominate the infiaton potential, we find = 0.976(1), 
where the uncertainty comes from the spread of reheat temperatures in that range. 
The cosmic string mass per unit length is greatly reduced by the presence of a Higgs 
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condensate at the string core, and is rendered independent of the inflaton couphng. 
Finally, thermal inflation sweeps away the gravitinos generated at the first stage of 
inflation, and any GUT-scale relics such as magnetic monopoles. 

There are other models which renormalisably couple F-term hybrid inflation to 
the MSSM. Fd hybrid inflation [14, 15] has the same field content as ours, but the 
MSSM has no U(l)' charges; and it requires a Fayet Iliopoulos term. Also potentially 
in the class is the B— L model of Refs. [16-18], although there is no explicit discussion 
of the coupling of the inflaton to the Higgs fields. In the model of Ref. [19] the 
waterfall fields are SU(2)j:j triplets. The authors identified a flat direction involving 
the Higgs, without pursuing its consequences. The original F-term inflation model [1] 
had a spontaneously broken global U(l) symmetry, and models based on coupling it to 
the MSSM have recently been explored in [20], again without the possibility of Higgs 
thermal inflation being noticed. The same field content can also produce a promising 
superconformal D-term inflation model [21]. 

Further afield, it is also possible to construct renormalisable models of inflation 
in the Next-to-Minimal Supersymmetric Standard Model using soft terms to generate 
the vacuum energy [22]. Inflation along a flat direction which mixes a singlet with an 
MSSM flat direction has also been investigated recently in Ref. [23]. In that work, a 
single stage of inflation was envisaged, and in order to supply a satisfactory spectral 
index, the coupling to the inflaton has to be non-renormalisable. 

The spectral index problem can also be solved with a non-minimal Kahler po- 
tential [24], or tuning the inflaton coupling to be small enough that the linear soft 
term dominates its potential [25]. In this paper we will restrict ourselves to the case 
where radiative corrections dominate the inflaton potential, and the Kahler potential 
is canonical. 

2 Coupling F-term inflation and the MSSM 

Our guiding principle is to couple minimal F-term hybrid inflation and the MSSM 
(which we take to include 3 families of right-handed neutrinos) in a renormalisable 
way, preserving all symmetries including supersymmetry (while allowing soft breaking 
terms in both sectors). Hence the superpotential will take the form 

W = Wi + Wa + Wx (2.1) 

where Wi is the standard linear F-term hybrid inflation superpotential of Eq. (1.1), 
Wa is the MSSM Yukawa superpotential 

Wa = H^QYuU + H^QYdD + H^LYeE + H^LY^N, (2.2) 

and Wx is the coupling between the inflaton sector and the MSSM superpotential, 
containing renormalisable terms only. We will assume that the U(l)' symmetry of the 
waterfall fields 

$ ^ $' = e^^*^$, $ ^ $' = e^^^$ (2.3) 
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U D 


Hi H2 N 






-Qe - Iql Qe + ^Ql 


-Qe -qi qE + Ql -2gL - Qe 



Table 1. Anomaly free U(l) charges for lepton doublet, singlet charges ql, qe respectively. 

is gauged. The inflaton S must be a gauge singlet, and so = — The mass scale 
M sets the inflation scale and the vevs of $ and $. Given that the inflation scale is of 
order 10^^ GeV, the waterfall fields must be SU(3)(8)SU(2)(g)U(l)y singlets_Note that 
Wi has a global U(l) R-symmetry, which forbids the terms 5^, 5*^ and $$. In order 
to preserve the flat potential for the inflaton, we must preserve this symmetry; we will 
discuss more of its implications in a moment. 

The form of Wx is now tightly constrained by symmetry and anomaly cancella- 
tion. Possible anomaly-free U(l) charge assignments for the MSSM fields are shown in 
Table 1. The SM gauged U(l)y is qi = —i,qE = 2. U(l)^__^ is qE = —qi = 1; in the 
absence of this would have U(l)^ and U(l)-gravitational anomalies. The diagonal 
subgroup of SU(2)r is qL = 0,qE = 1. Note that quite generally g^i = —qH2, so we will 
write qH2 = —qHi = qn- We will assume that the MSSM fields couple to a U(l)' distinct 
from U(l)y, i.e. that 2qL + qE 7^ 0, and moreover that in the AMSB case the values of qi 
and qE result in a solution to the AMSB tachyonic slepton problem [26]. For the result- 
ing sparticle spectra in this case, see Ref. [13]. (Note that if the U(l)' does not couple 
to MSSM fields, we are driven to Fd inflation [14, 15]). Three SU(3)(8)SU(2)(8)U(l)y 
singlets quadratic in the MSSM fields are available for Wx, namely H1H2, LH2 and 
A^A^ [27]. The U(l)' charge assignments, combined with the global R-symmetry, with 
superfield charges 

S = 2,L = E = N = U = D = Q = l,Hi = H2 = ^= ^ = 0, (2.4) 

now uniquely specify the coupling term as 

Wx = lX2NN<i> - \3SHiH2, (2.5) 

where we have set = ±(4^^, + 2qE) to permit the first term. All renormalisable B, 
L violating interactions and the A^A^ and LH2 mass terms are forbidden by the U(l)' 
gauge invariance, and the superpotential Eq. (2.1) contains all renormalisable terms 
consistent with U(l)' and the R-symmetry. Note in particular that the R-symmetry 
forbids the Higgs //-term H1H2. Moreover, the R-symmetry forbids the quartic super- 
potential terms QQQL and UUDE, which are allowed by the U(l) symmetry, and 
give rise to dimension 5 operators capable of causing proton decay [28, 29]. In fact the 
charges in Eq. (2.4) disallow B-violating operators in the superpotential of arbitrary 
dimension. 

Soft terms break the continuous R-symmetry to the usual R-parity. The lightest 
supersymmetric particle (LSP) is therefore stable. (From Eq. (2.4), the LSP is a scalar 
quark or lepton, or a gaugino, or a fermionic Higgs, S, $ or $.) 
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To summarise the assumptions which force us to this unique class of theories, we 
require a theory with : 

1. The field content of minimal F-term inflation and the MSSM. 

2. The symmetries of minimal F-term inflation and the MSSM. 

3. Renormalisable couplings only. 

4. An inflaton-sector U(l)' gauge symmetry which is coupled to the MSSM. 

Note that if $ and $ are gauged under a larger symmetry group, the coupling 
iViV$ is not allowed, unless they are triplets of SU(2)^ and {N,E) are doublets [19]. 

The parameters M, Ai, A3 are real and positive and A2 is a symmetric 3x3 matrix 
which we will take to be real and diagonal. The sign of the A3 term above is chosen 
because with our conventions, in the electroweak vacuum 

we have H1H2 — t- —\viV2- 

In the following we will denote the SU(3)®SU(2)(8)U(l)y gauge couplings by 5^3, 
g2 and gi, and the U(l)' gauge coupling by g'. The normalisation of the U(l)y gauge 
coupling corresponds to the usual SM convention, not that appropriate for SU(5) uni- 
fication. We will denote the soft parameters for the gaugino masses Ma, for a cubic 
interaction with Yukawa coupling A hx, and for a mass term 0*0 (where denotes a 
scalar field), m|. For the one mass term of the form 0^ in the MSSM {H1H2) we will 
use ml. 



3 The Higgs potential and its extrema 

In this section we explore the important extrema of the Higgs potential, and demon- 
strate that there is a 1-parameter family of supersymmetric ground states with non-zero 
vevs for 0, and hi^2 before supersymmetry-breaking is taken into account. We will 
assume that M, the scale of inflation and U(l)' symmetry-breaking, is much larger 
than the scale of supersymmetry-breaking. 

The existence of the one-parameter family (before thermal effects and soft terms 
are taken into account), is demonstrated as follows. The minimum of the scalar po- 
tential is determined by the requirement that both the F- and D-terms vanish. The 
vanishing of the D-terms ensures that |0| = |0|, \hi\ = I/12I and h\h2 = 0, while the 
vanishing of the F-term is assured by Ai00— A3/11/12 = M^. The minimum can therefore 
be parametrised by an SU(2) gauge transformation and angles x? defined by 

M 

(hi) ~ ia2(h2Y ^ (^^cosx,0), 

V A3 

— * M 

(0)^(0)^ sin xe*^. (3.1) 
V Ai 
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The cp angle can always be removed by a U(l)" gauge transformation (where the 
residual symmetry unbroken by the Higgs vevs alone is 11(1)^^^ x U(l) ), so the physical 
flat direction just maps out the interval < x ^ 1^/2. At the special point x = the 
U(l)" symmetry is restored, and at x = 7r/2 the SU(2) (8)U(l)y is restored. Away from 
these special points only U(l)em is unbroken. 

The degenerate minima have been noted before [19] in a model with gauge group 
SU(3) ® SU(2)^ ® SU(2)^ ® U(l)^_^. However, the important cosmological conse- 
quences which follow was first explored in Ref. [13]. 

Let us first consider the limiting cases where either hi^2 or 0, vanish. 



3.1 The (f),4>,s extremum (0-vacuum) 

In the 0, 0, s subspace (lower case fields denote the scalar component of the superfields) 
the scalar potential (including soft supersymmetry-breaking terms) is: 



V = A?(|0s|2 + |0sn + |Ai^ 



2|2 I 1^2/2 



-T,')\2 



+ mj|0p + m||0|^ + m^|s|^ + pM'^m^{s + s*) 
+ hxMs + c.c. 



(3.2) 



We will assume that the term linear in s is small enough not to be important for 
inflation (and quantify this smallness in Section 5). In AMSB there are arguments [30] 
to show that, without a quadratic term S*^ in the superpotential, the only RG invariant 
solution for p is p = 0. 

Let us establish the minimum in this subspace, under the assumption that ma <C 

M. We shall call this the 0- vacuum. With the notation (0) = v^j^^^ (0) = i^/ \/2 

and (s) = Vs/ -\/2, we find 



ml + IXiivl + v^)] + ^v^v^ + V2pM 



rris 

2 



From Eqs. (3.3), (3.4) we find 



'Z , I , \2 2 

ms + mj + X^v 



•^ml-vlml + {vl-vl)lXlv^ 
vl + vl 



0, (3.3) 
0, (3.4) 
0. (3.5) 

(3.6) 
(3.7) 



Then from Eqs. (3.6), (3.7), to leading order in an expansion in mzjM we have 



and from Eq. (3.5) that Vs is 0(m3). It follows from Eq. (3.7) that 



2 2 



9'M 



+ 0{mi/M^) 



(3.9) 



and from Eq. (3.5) that 



2 



V2XI V2Xi 
From now on we neglect p, assuming that 



0(m|/M). 



(3.10) 



\P\ 



< 



Aims 



(3.11) 



Substituting back from Eqs. (3.8), (3.10) into Eq. (3.2), we obtain to leading order 



Ai 



and from Eq. (3.10) a Higgs /x-term 



2\i 



(3.12) 



Ph 



2A? ' 



(3.13) 



naturally of the same order as the supersymmetry-breaking scale. 

The theory is approximately supersymmetric at the scale M, so the U(l) gauge 
boson, the Higgs boson, the gaugino and one combination of t/;^^ form a massive 

supermultiplet with mass m ~ g' ^Jv^^'+v^, while the remaining combination of (p and 
(f) and the other combination of tp^-^ form a massive chiral supermultiplet, with mass 

The large vev for generates inflation-scale masses for the triplet, thus natu- 
rally implementing the see-saw mechanism. 

3.2 The /ii 2, s extremum (/i-vacuum) 

. In the hi^2, s subspace, the scalar potential is 



2\2 



V = Xl{\h,s\' + IMP) + \X,h,h2 - M2|2 + Ig'^ql [\h,\^ - \h,\^) 

a 

+ m^J/iip + m^J/i2p + + pM'^mz^s + s*) 

+ hx.ihih2S + c.c. 



(3.14) 



detailed explanation of the symmetry-breaking is contained in Ref. [18]. 



-8- 



Note that we assume there is no hih2 mass term; its absence follows from the absence of 
the corresponding term in the superpotential (which is forbidden by the R-symmetry) 
when the source of supersymmetry breaking can be represented by a non-zero vev for 
a spurion (or conformal compensator) field. 

The structure is similar to Eq. (3.2), with the addition of SU(2) and U{1)y D- 
terms. Without loss of generality the SU(2) D-term vanishes with the choice hi = 
(wi/a/2, 0) and h2 = (0,W2/v^), and vi = V2- The values of the fields at the minimum 
(which we term the /i- vacuum) and the value of the potential at this extremum can then 
be recovered from the result of the previous section with the replacement Ai — A3), 
leading to a potential energy density 



3.3 Potential along the 0, </>, hi, h2 fiat direction 

As we outlined at the beginning of the section, the supersymmetric minima are parametrised 
by an angle defined in (3.1). Soft terms lift this degeneracy, and the leading terms 

in the effective potential for x can be found in an expansion in m\/M'^. After solving 

2 

for s, it is found that 

M2 f^A^sin^x + ZiAaCos^x) 2^m2 m? 2 \ . x 

V{x) ^ -— ^ X ■ 2 , , + Msin^x + -^cos^X , (3.16) 

2 Ai sm X + A3 cos"^ x \ Ai A3 / 

where we have defined 

I r h\,^ _ 2 2,2 -2 2,2 ^01 

Ai A3 

4 Supersymmetry- breaking and the true minimum 

In this section we investigate under which conditions the phenomenologically accept- 
able large-0 solution is the true minimum, in three popular supersymmetry-breaking 
scenarios. Hence we are looking for constraints on the soft supersymmetry-breaking 
parameters such that 

Vh-V, = M'\'^-'^-'^+'^\>Q, (4.1) 



™2 



V^"(V2) = ^ ^\2^-^-l\+mi-^-m 



> 0. (4.2) 



We will also check that the false vacuum at x = is a local maximum, from the 
sign of V"{0), which can be recovered from V"{7r/2) by the replacements 1^3 and 
-H- m|. A metastable false vacuum, as we will demonstrate in Section 7, would 
lead to the universe remaining trapped in an inflating phase. 
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We assume that the U(l)' symmetry is broken by a vev of order ~ M/ ^/Ai^ , 
and evaluate the soft terms at this scale, rather than running down to the electroweak 
scale. This is the appropriate renormalisation scale to investigate a potential with 
vevs of order w', whose important radiative corrections are from particles of mass of 
order g'v' and M . Note that in inflation models, with infiaton couplings Ai and A3 are 

generally small, and so the U(l)' gauge boson mass niA = g' ^Jv^^'+v^ is much greater 
than M, unless g' is also small. 

4.1 Anomaly-mediated supersymmetry-breaking 

With anomaly mediation, the soft breaking parameters take the generic renormalisation 
group invariant form 

Ma = ms^gjga, (4.3) 

hu,D,E,N = -m^PYu,D,E,N^ (4-4) 

{^yj = l^f^^l'^ + kY:6^v (4.5) 
ml = Km-A^h ~ TTT'3f3^^. (4.6) 

Here fi is the renormalisation scale, and ma is the gravitino mass; are the gauge 
/3-functions and 7 is the chiral supermultipfet anomalous dimension matrix. Yu^£),e,n 
are the 3x3 Yukawa matrices, fih is the superpotential Higgs /i-term, n and k are 
constants, and Y- are charges corresponding to the U(l) symmetry. 

In the MSSM, k is an arbitrary parameter, which in practice is fixed by minimising 
the Higgs potential at the electroweak scale. The parameter k is generated by the 
breaking of the U(l) symmetry at a large scale, and forms the basis of the solution 
to the tachyonic slepton problem within the framework of AMSB, as explained in 
[13], whence the name strictly anomaly-mediated supersymmetry-breaking (sAMSB) 
originates. 

The Higgs /i-term, fi^, is generated by the the vev of the inflation s, which in turn 
is triggered by the U(l)' symmetry-breaking. Hence the parameter k, and the equation 
for m|, are relevant only below the U(l)' symmetry-breaking scale v'. 

As a first approximation, we will assume that the g' terms dominate throughout, 
as qn and g$ are generally large, in which case the hx^ and trilinear soft terms are 
given from Eq. (4.4) as: 



hx, ~m3-i- 4g|/ , (4.7) 



2 I67r2 

^A3 ^ ma A_(4g|/), (4.8) 
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while the mass soft terms are 



d 

m=r ~ —7723 



2 _ _2 ^ 



2 2 



2 2 



mr ~ -rris 



The one loop g' /3-function is 



/^.' = (4-13) 



where 

Q = ncifql + 8g| + IGg^jgi) + 36ql + AOqEqi + 12g| 

= 76gi + 36g| + 88g£;gi (4.14) 

for riG = 3. Hence 

^ m| ^ -2m| j g|Q, (4.15) 

/ /2 N 2 

< ^ < ^ -2r«| (^^J g^Q- (4-16) 

Thus the difference in the energy densities between the two vacua is, in this approxi- 
mation, 

,/2\ 2 , 



Vh-V^-M 



2 f 2 



m9-\ HQql + 8ql ^Qqjj + 8qj, 



167r2 



Ai A3 



(4.17) 



The coefficient Q is in general large, and larger than both g| and qjj, so the condition 
for Vrp to be the true minimum may be written 



Ai ~ V?* 



(4.18) 



It is not hard to check from Eq. (4.2)) that under the same assumptions, the 0- vacuum 
is a minimum and the /i-vacuum is a maximum. Hence no further constraints on the 
parameters are generated. 

In the next section we will see that if A3 > Ai, then inflation ends with 0,0 
developing non-zero vevs, whereas if A3 < Ai it is (/ii,2) which become non-zero; this 
statement is independent of the nature of the soft breaking terms. Now is easy to show 

that < 1 unless 



9* 



- - < ^ < -1. (4.19) 



-II- 



However, the domain defined by Eq. (4.19) does not permit a satisfactory electroweak 
vacuum in tlie AMSB case [26]. For example, for the specific choice ql = 0, which can 
lead to an acceptable electro-weak vacuum [13], the condition Vh > becomes (from 
Eq. (4.17)) 

or Ai < 4A3 from the approximation Eq. (4.18). 

We see, therefore, that there will generally be a domain 

such that the universe exits to the false high Higgs vev /i- vacuum, evolving subsequently 
to the true vacuum as we shall describe later. 

In the Appendix we include a more accurate computation of the vacuum energy 
difference, taking into account the SM gauge couplings and the top Yukawa coupling. 

4.2 Gauge- mediated supersymmetry-breaking 

In the GMSB framework (see e.g. [31]), supersymmetry-breaking is communicated by a 
set of messenger fields C which have SM gauge charges in a vector-like representation, 
which should be complete GUT multiplets if gauge unification is to be preserved. The 
messenger fields are supposed to have a large mass, given by the vev of the scalar 
component of a chiral superfield X, which also has a non-zero F-term Fx-, the source 
of the supersymmetry breaking. Although there are many possibly choices for the field 
representations of the messenger fields, we can adapt the simple model described in 
[31] to study our model. 

We introduce the following superpotential for the extra fields 

Wg^ = X^SCC + X5XCC, (4.22) 

assuming that some extra dynamics at a higher scale gives both the scalar component 
of X and Fx a vev. We will assume that (X) ^ M. Radiative corrections from the 
messenger particles then induce masses for the gauginos at one loop, 

= ji^A,, (4.23) 

where Ag = N^i{Fx)/Mx, Mx = X5{X), and N^i is the messenger index, equal to 
twice the sum of the Dynkin indices of the messenger fields. Scalars acquire masses 
from 2-loop corrections of 



-.2 \ 2 



where A^ = N^i{{Fx) / MxY , Ca{i) is the quadratic Casimir associated with the ath 
gauge group for the ith scalar, and the sum over a includes the four gauge couplings 
91^2., 9' ■ 
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Trilinear terms are also induced at 2 loops, and so are of order Ag(aa/47r)^. They 
are small compared with the gaugino masses, and it is a reasonable approximation to 
take them to vanish at the messenger scale Mx- We assume that A^^^ are of the correct 
order of magnitude for supersymmetry-breaking. 

We thus have 



Ai 



m- 



2A 



9 



i2 



As 



3 



9l 



4 V 167r2 



9i_ 
4 V 167r2 



,'2 



167r2 



(4.25) 



Thus the difference between the vacuum energies is 



2A2M2 
(167r2)2 



As 
3 



ml + m| 

aT 



9t 



-A 



2ql9" 



1 

A^ 



2 



Ai 



(4.26) 



2"^ 2^ 

so that, if we assume dominance of the g' terms, the condition that < Vh becomes 

2 



As< 
Ai ~ 



qH_ 



(4.27) 



This is precisely the opposite condition to that in AMSB, Eq. (4.18). As in AMSB, the 
condition that V"^ < \4 is sufficient to ensure that is a minimum and Vh a maximum. 

Now in GMSB, we do not have the constraint on the domain {qliCIe) that we 
described in the AMSB case. Inflation will end in the Higgs phase unless 



qH_ 



> 1 and 1 < ^ < 
Ai 



qH_ 



(4.28) 



in which case it ends directly in the true 0-vacuum. 

4.3 Constrained minimal supersymmetric standard model 

At the high scale we will have the CMSSM pattern of soft breaking parameters. 



^1 


<p 




hxs 


Ai 


As 



A 



m 



h2 



01 



and hence 



Vh-V^ = M\2ml-A^/2) 



1 

A^ 



1 

A^ 



(4.29) 



(4.30) 
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Hence if A3 < Ai (so that inflation ends in the /i-vacuum) then for Vh > we require 

2ml > ^V2- (4.31) 

It is easy to check from Eq. (4.2) that this is again a sufficient condition that V"{7r/2) 
be positive. On the other hand, there is then a range 

for which the /i-vacuum is also a local minimum. We will see that this scenario is not 
consistent with a graceful exit from Higgs thermal inflation, and hence for a cosmolog- 
ically acceptable potential, we must demand 

> (4.33) 

5 Inflation and reheating 
5.1 F-term inflation 

We assume that the vevs of MSSM fields apart from the Higgs are negligible, in which 
case the relevant tree potential is 

Vt,ee = |Al00 - Xshh2 - MY + [A?(|0|2 + |0|2) + Xl(\h,\^ + j/^^j^)] 

^g'"^ - 0*0) + <lH{h\hi - h\h2^ 



+ hi + hl^''h2)' + hMhi - hlh 



2/ 



+ V^soft. (5.1) 

The soft terms in Vgoft are those appearing in Eqs. (3.2), (3.14), and are all suppressed 
by at least one power of ms. The most important soft term for inflation is one linear 
in s, the effect of which we assume is small compared with the radiative correction. 
We will see in Eq. (5.7) that this implies tuning below 0(1) only if the couplings Ai^s 
are very small. We also assume that the higher order terms in the Kahler potential do 
not contribute significantly. 

At large s, and with all other fields vanishing, the potential is approximately 

V = M^ + AVu (5.2) 

where AVi represents the one-loop corrections, which dominate the soft terms. As S 
is coupled only to $, $ and Hi 2, the contribution to the one-loop scalar potential is 

[32] 



32n 



2 



(Xjs^ + A,M=)= In ('llif-tj^iM!^ + - A,M=)Mn ('MfizilM! 



X2e2\ /\2„2 



-2Al.Mn(-^J-4A^.Mn^^, 



(5.3) 
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For large s (meaning Ai^ss^ ^ M^) the potential can be written as 



V{s) - 



1 + a In ■ 

si 



(5.4) 



where an 0(a) correction to has been dropped, and 



2 



a = A_, A=^A? + 2Ai, si = M''/\. (5.5) 

We will neglect supergravity contributions in the potential, which will require a small 
coupling c of the quartic term c|s|^/mp in the Kahler potential, and impose a constraint 



A < 0.06. (5.6) 

There are also potentially important contributions from the linear soft term pM^ma s + 
c.c. These are negligible provided 

P « (5-7) 
167r^ ms 

2 

We will shortly see that Sc ~ 10^® GeV, so assuming ms ~ 10^ GeV, a soft term with 
p ~ 1 is negligible provided 

10-^ < A. (5.8) 

Henceforth we will assume that the Kahler potential is canonical and that A is in the 
range given by Eqs. (5.6), (5.8). We note, however, that interesting consequences for 
the spectral index flow from a non-canonical Kahler potential [24] and from couplings 
small enough for the soft term to contribute [25]. 

5.2 Perturbation amplitudes 

The scalar and tensor power spectra Vg-, Vt and the scalar spectral index Ug generated 
on a scale k equal to the co-moving Hubble scale aH at A^^ e-foldings before the end 
of inflation are given by the standard formulae (see e.g. [34]), 

f.W = Trfr— I — I =^|— V. (5-9) 









3 






' 8 ( 


Sc \ 








nip J 



r^s^il-^]. (5.11) 



The WMAP7 best-fit values for Vsiko) and at a pivot scale k = ko = 0.002 /iMpc ^ 
in the standard ACDM model are [35] 

Vs{ko) = (2.43 ± 0.11) X 10-^ Us = 0.963 ± 0.012(68%CL). (5.12) 
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From this data we infer that 

— ^ 2.9 X 10-3 fl^y , N,, = 2lt\\ (5.13) 

showing approximately a 2a discrepancy with the standard Hot Big Bang result N^^ ~ 
58 + ln(Trh/10^^ GeV) (assuming only MSSM degrees of freedom at T^h)- We will see 
shortly that the reheat temperature lies in a range around 10^^ GeV, and in Section 7 
that there are Ng ~ 15 e-foldings of thermal inflation at a lower scale. Therefore one can 
estimate A^pti — 42(1) e-foldings of F-term inflation while the pivot scale ko is outside 
the horizon, where the uncertainty comes from the range of reheat temperatures, given 
in Eq. (5.19). The scalar spectral index is thereby reduced to 

~ (l - ~ 0.976(1). (5.14) 

V ^Fti / 

Lower values of the spectral index are possible if A drops below the limit (5.8) and the 
linear soft term comes into play [25]. 

5.3 End of inflation and reheating 

F-term inflation ends when one set of scalar fields becomes unstable. If A3 > Ai, the 
(j), (f) pair become unstable first, and inflation ends at the critical value s^^ = M^/Ai. 
The fields 0, gain vevs and the universe makes a transition to the U(l) -broken 
phase described by Eq. (3.3)-Eq. (3.5). On the other hand, if A3 < Ai, the Higgs fields 
become unstable first, the critical value of s is 5^3 = MVA3, and the universe makes a 
transition to a phase where hi and /12 develop vevs of order the unification scale rather 
than 0, 0. In this phase the SU{2)l symmetry is broken. 

At first sight, this would appear to rule out the model with A3 < Ai. However, 
provided the correct (small Higgs vev) vacuum has the lowest energy density at zero 
temperature, the universe can seek the true vacuum when thermal corrections become 
sub-dominant. We will establish in Section 7 that the evolution to the true ground 
state proceeds by a period of inflation. 

Assuming that A3 < Ai, inflation exits to the /i- vacuum, with symmetry-breaking 

SU(2)®U(1)^®U(1)' ^ U(1)^^®U(1)". (5.15) 

Here, U(l)" is generated by the linear combination of hyper charge and U(l)' generators 
which leaves the Higgses invariant: 

Y" = Y' -{qL + qE)Y. (5.16) 

There are still two Abelian symmetries, and SU(2) is completely broken with no discrete 
subgroup. Hence cosmic strings are not formed at this transition. 

We expect reheating to be very rapid [36-41], as the period of oscillation of the 
fields is of order M~^, which is much less than a Hubble time, and the couplings of 
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the Higgs field are not all small. Hence the universe regains a relativistic equation of 
state almost immediately, and thermalises at a temperature Trhi given by 



» VM.f^VVA,. (5.17) 



where (?rhi is the effective number of relativistic degrees of freedom at temperature Trhi. 
From (5.13), and taking ^frhi = 915/4 (a slight overestimate), we find 

T,hi - 2.2VA X 10^^ GeV. (5.18) 

Hence the range of reheat temperatures corresponding to the range of couplings defined 
by Eq. (5.6) and Eq. (5.8) is 

0.7 X 10^^ < Trhi/GeV < 5 x 10^^ (5.19) 

Finally, we note that large vevs of other fields along supersymmetric flat directions 
can lead to blocking of particle production during reheating [27]. On the other hand, 
radiative corrections during inflation generically generate masses of order y^H'^ [42], 
where y is a combination of Yukawa couplings, and so we expect that other vevs besides 
that of the infiaton will be generally small. We leave a detailed examination of the fiat 
directions for another work, assuming for now that any fiat directions which do not 
have y of order 1 are small. 

5.4 High temperature ground state 

As the universe reheats, it will seek a minimum of the finite temperature effective 
potential, or equivalently the free energy density. To discuss the free energy, it is 
convenient to define a dimensionful field X = f +x, with v+ = V2V(j, = 2M / The 
free energy density can then be expressed as 

/(X,T) = -^^7efr(X,T)T^ (5.20) 

where g^siX^ T) is the effective number of relativistic degrees of freedom at temperature 
T. At weak coupling, gcsi^X, T) can be calculated in the high-temperature expansion 
for a particle of mass m '^T [43], 

90 

(7eff(^,T)~Co-Ci — — , (5.21) 

where there are contributions to cq of 1, | and to Ci of ^, ^ for bosons and fermions 
respectively. For particles with m > T, g^fi is exponentially suppressed. 

We can see that X = is a local minimum for temperatures ma ^ T < M, 
because away from that point the U(l)" gauge boson develops a mass proportional to 
(0), and so ges decreases. For similar reasons X^ = v+7t/2 is also a local minimum: 
away from that point the MSSM particles develop masses and again reduce ges- 
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In fact, by counting relativistic degrees of freedom at temperatures ms <^T < M 
one finds that is the global minimum at high temperature. In the /z-vacuum the 
relativistic species are the $, $ chiral multiplets and the U(l)" gauge multiplet. In the 
0-vacuum, the particles of the MSSM are all light relative to T. Hence 

/(0,T)^-y^n (5.22) 

f{X,,T)c^- — -T\ (5.23) 

The minima of the free energy density are separated by a free energy barrier of height 
~ T^. The transition rate can be calculated in the standard way [44] by calculating 
the free energy of the critical bubble Ec. The transition rate per unit volume is then 

|)'e.p(-f). (6.24, 

The critical bubble is a solution to the equation 

X" + -X' + V^siX) = 0, (5.25) 
r 

where r is the radial distance from the bubble centre, and we have neglected 0(1) 
complications in the kinetic term from the non-linear field transformation. An order- 
of-magnitude estimate can be given, recognising that X has to change by an amount 
AX ~ f+ from the inside to the outside of the bubble, while negotiating a local free 
energy bump of order AV^ ~ T^. Neglecting the damping term, one can translate the 
equation into a harmonic oscillator problem, finding that the critical bubble radius is 
approximately 

AX 



(5.26) 



and so the critical bubble energy 



~ AV>3 ~ -4£L ~ ^. (5.27) 

The universe will stay in the wrong ground state if the transition rate per unit volume 
is significantly below the Hubble rate per Hubble volume, or F < H'^. Hence the reheat 
temperature T^hi should be parametrically 

Ta < — r- (5.28) 

[ln(mp/u+)] 3 

Recalling that Trh — M and v+ = 2M/ \/Xi, we see that if inflation exits to the h- 
vacuum it is likely that the universe continues to evolve with large (inflation-scale) 
Higgs vevs, provided Ai ^ 1. 
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6 Review of gravitino constraints 



There are strong constraints on the gravitino mass and hfetime from cosmology [9- 
12]. If the gravitinos are unstable, they can conflict with Big Bang Nucleosynthesis 
(BBN) by photodissociating light elements, or they can decay directly into the LSP, 
which in turn produces a limit from the known density of dark matter in the standard 
cosmological model. The gravitino may also be the LSP, in which case the dark matter 
constraint applies directly. 

Gravitinos are produced by collisions of high-energy particles in the thermal bath, 
principally gluons and gluinos, with an abundance of approximately [45] 

, \ / T 

-th _ , , / o /I I 1 /I 9 \ in-13 / 



I 2,4 + .4^ )x 10- ^^ ). (6,i; 



2 



where Mg is the gaugino mass at the GUT scale. We include an 0(1) factor Uq to take 
into account the theoretical uncertainties [46-48], arising from the strong dynamics of 
the coloured plasma. 

BBN constraints [45] are not easily summarised, but are much tighter for lighter 
gravitinos which decay during or after nucleosynthesis, as relevant for the CMSSM. For 
gravitino masses less than about 0(10) TeV, the reheat temperature is bounded above 
by ^rh ^ (0.2 — 1) X 10^ GeV. For higher gravitino masses, the dark matter density 
provides a bound, and so it is appropriate use Eq. (6.1) in the limit M?/m| — )■ 0.. 

^ 2 

Given that the LSP density parameter arising from a particular relic abundance 
in the MSSM is 

the LSP density parameter from (high mass) thermally produced gravitinos can be 
found as 

This must be less than or equal to the dark matter abundance inferred from the CMB 

[35] 

fidm/i^ ^ 0.11. (6.4) 

The presence of cosmic strings in our model, although affecting the CMB power spec- 
trum, does not significantly affect this inferred value [49]. 

In our model, we will see that the gravitinos generated by the first stage of reheat- 
ing are diluted by a period of thermal inflation. The constraint therefore applies to 
reheating after thermal inflation. We will also see that the second reheat temperature 
is about 10^ GeV, and so we can only tolerate unstable gravitinos of mass greater than 
about 10 TeV in order not to spoil BBN. This is natural in AMSB, problematic in 
GMSB, while the CMSSM keeps ma as a separate parameter. 

2 

There are also non-thermal production mechanisms from coherent oscillations of 
the inflaton [50, 51] and from ordinary perturbative decay [52], whose rates depend on 
the inflaton mass and vev. We will see in the next section that the relevant inflaton 
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mass and vev will be those of the Higgs. However the BBN constraints mean that the 
gravitino, when it is not the LSP, must be much more massive than the Higgs and so 
cannot be produced by direct decays. Hence only thermal production is relevant. 



7 Higgs thermal inflation and gravitinos 

As the temperature falls, the energy density difference between the vacua becomes 
comparable to thermal energy density, and the universe can seek its true ground state, 
which is X = ^/2, the 0- vacuum. 

At zero temperature we can write the difference in energy density between the 
/i-vacuum and the 0-vacuum as (see Eqs. (4.17), (A. 13), (4.26) and (4.30)) 



AV^ff - vim" 



sb 



(7.1) 



where we recall that v\ = AM^/Xi, and we have defined an effective SUSY-breaking 
scale rrtsb- In the supersymmetry-breaking scenarios under consideration 



ml 



A3 



9l 



8 \ 167r2 



A period of thermal inflation [53] starts at 



(AMSB), 
(GMSB), 
(CMSSM). 



(7.2) 



Ti 



30 2 2^^ 



(7.3) 



where g\ is the effective number of degrees of freedom at temperature T^. The CMB 
normalisation (5.13) for e-foldings of standard hybrid inflation gives {v+/m-p) ^ 
5 X 10~^(40/A^) 4 . Using the number of degrees of freedom for a ^iX)^^ theory 
with two light chiral multiplets $ and $, gi = 15, we have (on dropping the unimportant 
dependence on A^) 



Ti ~ 2.2 X 10" 



9 / rrish \ 2 



GeV. 



(7.4) 



.1 TeV/ 

The /i-vacuum must be a local maximum at zero temperature, i.e. the soft mass terms 
must be negative. If the /i-vacuum were a local minimum, one can 



P + mil 



estimate that the tunnelling rate per Hubble time per Hubble volume [54] would be 



r 

IF 



M8 



(7.5) 



where Se is the action of the Euclidean tunnelling solution. This ratio must be of order 
unity for the universe not to remain trapped in the false vacuum [55], and since the 
prefactor is much less than unity, we see that we cannot allow a metastable /i-vacuum 
for a graceful exit from thermal inflation. 
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Thermal inflation continues until the quadratic term in the thermal potential 
gf'^T^d^P + becomes the same size as the negative soft mass terms. Near the 
false vacuum, the high temperature effective potential for the field X breaking the 
U(l)" symmetry can be written [44, 56] 

Kfr(X) ^ i7(T2 - Ti)X' - \5TX^ + \\xX\ (7.6) 

where 7, 5 and Ax are dimensionless constants, and To ~ \m^\/g' . The cubic term 
arises from the gauge boson, and the transition is first order provided Ax < e^, where 
e is the effective U(l)" gauge coupling [44, 56]. 

Hence the transition which ends thermal inflation takes place at Te ~ TTisb, and 
the number of e-foldings of thermal inflation is 

^^-iln(2:)-15-ln(fl-), (7.7) 

Thus gravitinos will be diluted to unobservably low densities, as will any baryon number 
generated prior to thermal inflation, and any other dangerous GUT-scale relics such 
as monopoles. 

After thermal inflation ends, there is another period of reheating as the energy 
of the modulus X is converted to particles. Around the true vacuum, the X is mostly 
Higgs, and so its large amplitude oscillations will be quickly converted into the parti- 
cles of the MSSM. The natural oscillation frequency around x = is of order ma, 
while the Hubble rate is of order maM/mp. Hence in much less than an expansion 
time, the vacuum energy will be efficiently converted into thermal energy. The reheat 
temperature following thermal inflation is thus 



1 1 



where grh2 is the effective number of relativistic degrees of freedom at Trh2, given its 
MSSM value grh2 = 915/4. This second reheating regenerates the gravitinos, and we 
may apply the gravitino density formula Eq. (6.3), finding 

1 

ili^sph ~ 6 X 10 — (7.9) 

^ ^100 GeV VI TeV/ ^ ' 

We can convert the relic density into a constraint on the effective SUSY-breaking scale 
rrigb, requiring that the LSP density is less than or equal to the observed dark matter 
abundance, fldmh^ — 0.11. 

(J 

The parameter is directly related to physical observables differently in the different 
SUSY-breaking schemes, for which we can derive constraints. 
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7.1 Gravitino constraint in AMSB 

Using Eq. (7.10) and the expression for in Eq. (7.2), we find 

5 X lO'' / mrr,T.\ -2 



2 



Hence AMSB-based models requires a high gravitino mass in order to saturate the 
bound and generate the dark matter. 

We can be a bit more precise if we use use the phenomenological relations derived 
in [13]. Firstly, in order to fit fih we have (using Eqs. (3.13), (4.7)) 

dg" - ^, (7.12) 

while we can use a phenomenological formula for the LSP mass 

rriLSP ~ 3.3 x 10"^m|. (7.13) 

Hence 

We also have a constraint (Eq. (4.21)) on A3/A1 from requiring the exit to a false h- 
vacuum. Hence in order for the LSP in this model to comprise all the dark matter, we 
have 

1 1 



,^lVQ(lll ~ 360 TeV ^\ujI^ 



(7.15) 



For example, taking = as in [13], we find that ma is independent of and in the 
range 

130 TeV < ms < 250 TeV, (7.16) 

where we recall from the discussion around Eq. (6.1) that uJq is 0(1). It was noted 
in [13] that a Higgs of mass 125 GeV demands a gravitino mass of about 140 TeV in 
sAMSB, which is compatible with an LSP produced by gravitino decays being the dark 
matter. 

7.2 Gravitino constraint in GMSB 

In the GMSB framework the LSP is usually the gravitino, whose mass is given by 

m-A = —= A„, (7.17) 

where k < 1 parametrises the fraction of the total F-term contained in the messen- 
ger sector, and we recall that Mx is the messenger scale. It is more convenient to 
phrase the dark matter constraint in terms of the larger electroweak gaugino mass M2, 



-22- 



which dominates in the equation for the SUSY-breaking scale Ag, as in Eq. (4.23), and 
therefore Eq. (7.2) can be rewritten 



2 3 Ai M| 
" 8 A;iV~ 



(7.18) 



Hence 



^)%5^;i4(^)-. 

The bound can be saturated for a TeV-scale gravitino with TeV-scale gaugino masses 
without special tuning of the ratio A3/A1. A lighter gravitino forces the gaugino mass 
upwards. 

There is a separate constraint from decays of the NLSP (which is generally a 
neutralino for unless Mx is small), which may interfere with Big Bang Nucleosynthesis 
(see e.g. [31]). A careful analysis of the nucleosynthesis constraints [57] shows that a 
messenger mass of up to about 10^^ GeV is allowed, before hadronic jets injected after 
10'^ s results in the overproduction of ^Li. The combination of the dark matter and 
BBN constraints Mx < 10^^ GeV may be written 

(j^)'<0.5v/]^y5^. (7.20) 
7.3 Gravitino constraint in the CMSSM 

In the CMSSM, the gravitino bound Eq. (7.10) can be expressed in terms of the soft 
scalar masses mo and the trilinear parameter A from Eq. (7.2), as 

K - 5 X (I - i-O^y TeV. (7.21) 

This is a very weak bound, unless the ratio A3/A1 is very small: hence there is generi- 
cally a very low density of LSP dark matter generated by decays of gravitinos. Instead, 
the CMSSM can generate an acceptable dark matter density through the standard 
freeze-out scenario [58] (see [59-61] for recent analyses of the the CMSSM parameter 
space in the light of recent Higgs results). This requires that the gravitino mass is larger 
than about 10 TeV to avoid BBN constraints [45]. We conclude that the Higgs thermal 
inflation solution generally has no effect on the gravitino problem in the CMSSM, be- 
yond determining the reheat temperature and hence the standard BBN-induced lower 
bound gravitino mass. 

8 Conclusions 

In this paper we have shown how models which couple F-term hybrid inflation with 
the MSSM without extra ingredients naturally realise a period of thermal inflation, 
with a reheat temperature of around 10^ GeV, while generating the Higgs /x-term. The 
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inflation is driven by the relaxation of the Higgs fields to zero in a potential generated 
by the Higgs and waterfall field soft terms. This second period of intermediate scale 
inflation, which we have called Higgs thermal inflation, has a number of beneficial 
effects. It solves the gravitino overabundance problem of supersymmetric cosmology, 
while still maintaining the possibility of leptogenesis. It reduces the cosmic string mass 
per unit length so that CMB bounds are satisfied, and renders it independent of the 
inflaton couplings. Hence the scalar spectral index is not driven to unity in the effort 
to make the strings light, from which the tight constraints on standard F-term hybrid 
inflation are generated [8]. The period of thermal inflation means a reduced number 
of e-foldings of F-term inflation are required, and the scalar spectral index is reduced: 
in the range of inflaton couplings where the inflaton potential is dominated by the 
radiative corrections, we find Us — 0.976(1). 

The MHISSM is the simplest, attractive, formulation of this scenario, and gener- 
ates right-handed neutrino masses as well as the Higgs /i-term. We found constraints on 
the couplings and soft terms in order for the scenario to work: i.e. for F-term inflation to 
exit towards a vacuum with inflation-scale vevs for the Higgs field, and for that vacuum 
to be unstable. We investigated the implications of these constraints in three popular 
supersymmetry-breaking scenarios: AMSB (where the model coincides with strictly 
anomaly-mediated supersymmetry-breaking [13]), GMSB, and the CMSSM. We found 
constraints on the ratio of the inflaton couplings in AMSB (Eq. (4.18)) and GMSB 
(Eq. (4.27)), and that in the CMSSM the soft scalar mass must be greater than the 
half the magnitude of the soft trilinear term, multiplied by the square root of the ratio 
of the inflaton couplings (Eq. (4.33)). 

In AMSB, the gravitino problem becomes the gravitino solution: the observed 
dark matter density can be generated by the decays of gravitinos which are produced 
thermally following Higgs thermal inflation. In GMSB, the gravitino is the LSP, and a 
weak upper bound on its mass of about 1 TeV follows from the combined requirement 
that it supply the dark matter without NLSP decays spoiling nucleosynthesis. In the 
CMSSM, the density of thermally-produced gravitinos is generally sub-dominant, and 
the standard freeze-out scenario must do the work of making neutralino dark matter. 
However, the gravitinos must decay early enough not to spoil nucleosynthesis, meaning 
that the gravitino mass must be 0(10) TeV or greater. 

A reheat temperature of 10^ GeV is broadly consistent with thermal leptogenesis, 
provided at least one right-handed neutrino is light enough to be thermally produced. 
We leave the details for a future publication. 

The MHISSM predicts the formation cosmic strings, with dimensionless mass per 
unit length estimated as G/ig — 10"'^ [13]. While satisfying current CMB bounds, there 
are tight bounds on the GeV-scale cosmic 7-ray spectrum [62], so strings should have 
a very small branching fraction into 7. Strings may instead decay into gravitational 
waves, but there are also increasingly strict bounds on the stochastic gravitational 
wave background from pulsar timing [63, 64]. Should the bounds ultimately fall below 
the predicted value of G/ig, this will rule out the MHISSM but not the Higgs ther- 
mal inflation scenario in general, which remains a possibility whenever the inflaton is 
coupled to a set of waterfall fields which include the Higgs. 
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A sAMSB soft parameters and inflaton coupling constraints 

In this appendix we give formulae for a more accurate calculation of the sAMSB soft 
parameters in Section 4.1, relevant for the calculation of the vacuum energies and hence 
the constraint on A3/A1. 

The /lAi and h\.^ trilinear soft terms are determined in accordance with Eq. (4.4): 

hx, = -ms A_ (^sxl + ^ Tr + 2A^ - 4qlg"^ , (A.l) 
A, 



1x3 = -ms 



^— (Tr YeYI + 3 Tr Fzpi^z^ + 3 Tr YuVi + A? + 4 Tr A^ 



2 167r2 

-:igl-gl-Aqlg'^), (A.2) 
while the mass soft terms are 
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fA^ + Tr YeY^ + 3TtYbY^- 2g'\l - hf - M) , (A.5) 



1 d 
2n'^ ^ dfi 

IsL^^il + Trr;v4 + 3Trn.rJ - 2g"ql - \gl - \gi) . (A.6) 

In the following, we include the SM gauge couplings and the top Yukawa coupling, 
which we shall denote yt- We also retain the neutrino Yukawas Y^r, since their magni- 
tude is model dependent. We will assume that the value of tan/3 is such that it is a 
good approximation to neglect all the other Yukawas, and we will also neglect Ai and 
A3. 
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The relevant soft breaking parameters are then 
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We have assumed above for simphcity that, hke A2, F/^r is real and diagonal. 

The difference in the energies between the large-Higgs and large-0 vacua can be 
written 



(167r2)2Ai 



mg 

2 



2 ^2 - V^) ^ AQq^ + A„2 + 1 [Aqlg" + A,J 



Ai 



where 



A^2 



-i Tr A^ (1 Tr A2 + 4 Tr Y^yI - 2(g| + 2g^)(7'2) - Tr A 



-iTrAi 

-4 I Q„4 



(A. 13) 



(A. 14) 
(A.15) 



/2 



Ugt + 3^^ - 6 Tr Fj^ - 2 Tr F^^A 

-2 Tr F^ {Tr F^ + 3^^ - 3gl - gf _ 2(g^ + + g|,)(7 

- {(^yl - - 3^2 - fal - 2(gJ + + q^g" 
A,3 = gl + 3gl-3\yt\^-TTYj 



r2 

N 



(A.16) 
(A.17) 



- 26 - 



The condition on the couphngs deriving from the vacuum energies can therefore be 
written 

^ 1 + iqj,/Q) (/3|,A„. + 2 (1 + i/3^,A,3)') ^ 
^1 1 + (ql/Q) (/3|A^. + 2 (1 + i/3^A,j') ' 

where 

Taking the values of the SM couphngs at the U(l)' breaking scale to be the values at 
gauge coupling unification, we find using the renormalisation group analysis of Ref. [13] 

gi - 0.55, g2 ~ 0.71, ^3 ^ 0.70; yt ~ 0.51 (A.20) 

(where we have taken tan/3 = 16). Hence, at this level of accuracy, 

A,„2 ~ 3.5-6TrF^-2TrFj^A2 

-2 TrF^ {TrF^ - 1.0 - 2(g^ + ql + 

+ 1.6(gJ + g,^c + g|)/ (A.21) 

A;,3 ~ 1.0 -TrF^. (A.22) 

We can derive successive approximations. Firstly, neglecting terms of order qjj/Q and 
q^/Qi which is a good approximation given Eq. (4.14), we have 



^ > 4- (A.23) 
Ai 



Secondly, we can neglect terms of order I3h,4i and higher (which is not necessarily a 
good approximation), to obtain 

A3 g|l + 2(g|/g) 

Ai ql l + 2(g|/Q)- ^ ■ ^ 

In the case qi = 0, we obtain Eq. (4.20). 

We can also expand in powers of f3H,ii> while still neglecting terms of order Y^, and 
bearing in mind that an acceptable electroweak vacuum requires Tr A| ~ 4/3// [13]), we 
find to second order 

A3 qjj 1 + ^jqh/Q) (1 + [1-4 + OAjql + q%)/qjj]PH + 2.0/3|,) 

Ai^^l l + 2(g^/g)(l + lTrAi[l + 2gyg2]_i^^TrAi) ' ^ 
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